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ABSTRACT 


TJrls  repoffts  ^survey-y^^of  the  more  common  types  of  structural 
dynamics  problems  in  association  with  ballistic  missiles.  Special 
attention  is  given  to  the  effects  of  subsystem  interaction  on  dynamic 
response  of  the  missile  and  vice  versa.  A  few  technically  impor¬ 
tant  examples  are  treated  in  some  detail.  The  purpose  of  these 
examples  is  twofold;  first,  they  illustrate  various  techniques  used 
in  dynamics  analyses  and  second,  they  record  some  of  the  details 
of  mathematical  models  which  are  frequently  used  to  solve  missile 
dynamics  response  problems  Some  of  the  sections  of  this  report 
are  found  in  the  reference,afT^ut  are  repeated  here  to  provide  the 
reader  with  a  fairly  complete  pifcture  of  the  more  common  types  of 
missile  structural  dynamics  proMems. 
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1.  INTRODUCTION 


The  need  for  a  current  survey  and  reference  document  on  the 
more  common  miasile  structural  dynamic  problems  and  their  analyses 
resulted  in  the  writing  of  this  report.  Many  of  the  sections  of  the 
report  appear  in  previous  works,  as  referenced,  however  they  have 
been  revised  and  repeated  here  to  provide  an  overall  picture  of  the 
necessity  and  complexity  in  studying  missile  dynamic  response 
problems. 

The  assumption  of  a  rigid  structure  is  necessary  in  the  prelim¬ 
inary  design  phase  of  a  missile  in  order  to  allow  the  computation  of 
loads  on  which  further  design  work  can  be  based.  However,  once  a 
preliminary  structural  design  has  been  made  it  is  possible  to  compute 

the  dynamic  characteristics  of  the  missile.  Utilizing  these  dynamic 

* 

characteristics,  dynamic  loads,  a  revision  to  the  initial  load  estimates, 
can  be  determined. 

In  this  report  the  problem  of  determining  dynamic  loads  will  be 
discussed.  Special  attention  will  be  given  to  the  effects  of  subsystem 
interaction  on  dynamic  loads  and  vice  versa.  A  few  technically  impor¬ 
tant  examples  will  be  treated  in  some  detail.  ^ 

For  purposes  of  the  discussion  which  follows,  a  complete  missile 
system  is  envisioned  as  being  composed  of  the  following  subsystems; 

(A)  airframe  (structure)  including  propellants,  (B)  propulsion  system, 

(C)  control  system,  (D)  payload  (re-entry  vehicle,  satellite,  etc.  ), 

(E)  launch  facilities.  Each  of  these  subsystems  has  its  own  peculiar 
dynamic  characteristics  which,  from  the  dynamic  loads  point  of  view. 


Throughout  this  report  the  term  "dynamic  load"  is  used  to  designate 
those  stresses,  bending  moments,  shear  forces,  etc.,  which  are 
generated  by  a  time-dependent  environment.  The  stresses,  etc.,  are, 
of  course,  functions  of  time. 
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aro  not  exceptionally  interesting.  Hc'*-ever,  once  those  svibsystema 
are  interconnected,  the  picture  moL-r  change  completely.  For  example, 
a  propulsion  system  generates  a  thc-csi  which  builds  from  zero  to  a 
more  or  less  steady  valvie  and  there  liier  decays  to  coro.  The  time- 
dependent  character  of  the  thrust  p— ctinces  no  startling  behavior  when 
acting  on  a  rigid  airframe;  the  sifc.a:Cion  can  be  dramatically  different 
when  the  dynamic  behavior  of  the  airframe  is  considered.  These  typos 
of  dynamic  interactions  and  their  effect  on  dynamic  loads  are  the  cen¬ 
tral  theme  of  what  follows.  Some  c£  ihe  dynamical  environments  which 
are  considered  to  effect  the  missiler  response  are  ilhjstrated  in  Figure  1. 
The  effects  of  dynamical  environrnerrs  shown  in  Figure  1  which  are  not 
considered  in  this  report  such  as  uiraceady  aerodynamics  (bufleting)  and 
microseisms  are  given  in  References  1  to  3. 


Z.  IN- K  LIGHT  DYNAMIC  I, GADS 


In  this  section  equations  for  the  lateral  motion  of  a  liquid  pro¬ 
pellant  missile  will  be  develuped.  These  ec|uatit)iis ,  which  include  acM'D- 
dynamic  forces,  control  uysleiii  dynamics,  airframe  dynamics,  and 
propellant  dynamics,  can  be  used  to  compute  dynamic  loads  for  a  vari¬ 
ety  of  inputs.  In  addition,  they  can  also  be  used  to  study  control  system 
stability.  These  equations  are  evolved  in  References  4  through  12. 

2.  1  Airframe  Dynarnics  —  Lateral  Bendine 

Experimental  studies  have  established  the  fact  that  a  missile  air¬ 
frame  can  be  idealized  as  a  nonuniform  (in  stiffness  and  mass  distribu¬ 
tion)  beam  for  purposes  of  studying  its  lateral  bending  dynamics. 
Sometimes  it  may  be  necessary  to  represent  an  airframe  as  a  collec¬ 
tion  of  beams  rather  than  as  a  single  beam.  However,  such  a  represen¬ 
tation  poses  no  theoretical  difficulties,  although  computational  problems 
may  develop. 

It  has  been  established,  both  theoretically  and  experimentally, 
that  simple  (Bernoulli-Euler)  beam  theory  does  not  accurately  represent 
bending  dynamics  of  airframes  of  the  dimensions  encountered  in  ballistic 
missiles.  An  adequate  bending  theory  must  account  for  the  effects  of 
shear  deformation  and,  usually,  the  rotational  inertia  of  beam  elements. 
Such  a  theory,  known  as  the  Timoshenko  beam  theory,  is  available. 

The  Timoshenko  equations  are  derived  by  writing  the  equilibrium 
and  kinematic  equations  for  a  differential  element  of  the  beam  as  shown 
in  Figure  2,  The  beam  has  a  number  of  spring-mass  systems  attached 
at  distances  L  from  the  mass  center.  The  reasons  for  including  these 
systems  will  become  clear  later.  Equilibrium  of  the  lateral  forces  act¬ 
ing  on  the  element  requires  that  (see  References  5,  11  and  12) 


gives 


Figure  2.  Element  of  a  Timoshenko  Beam 


Writing  equilibrium  of  moments  and  dropping  second-order  terms 

=  .|M,v+p|2.  |^A(x.,j)=,^mjt.  (2! 


where  A^x  “  Dirac  delta  function  defined  by 

A^x  -  t  .\  =  0 


when  X  #  I  , 


and  ••vlu.TC 


u  =  total  dtiflection  of  bnam,  diio  to  both  bending  and  shear 
=  slope  of  beam  due  to  bending  alone 

=  dl fqjliifeinent  of  jth  attached  apriiig-inass  S'/sloni  relative 
■*  to  deflected  position  of  beam 

pA  =  the  mass  per  unit  length 

V  =  the  shear  force  on  a  cross  section 

pi  =  the  mass  moment  of  inertia  per  unit  length  of  the  beam 
cross  sections 

M  =  the  bending  moment  on  a  cross  section 
P  =  the  axial  load  in  the  beam  at  any  section 


a  =  absolute  acceleration  of  beam  system  in  negative 
x-direction 


=  number  of  slosh  masses  between  x  and  tip. 

Only  the  rotary  inertia  forces  caused  by  the  bending  slope  appear 
in  Equation  2.  The  shear  forces  produce  distortion  only  and  not  rota¬ 
tion  of  the  beam  elements.  The  moment  curvature  relationship  for  a 
Timoshenko  beam  is  the  same  as  for  a  simple  beam,  namely 


El 


(4) 


where  El  is  the  flexural  rigidity  of  ll\e  beam.  Because  of  the  inclusion 
of  shear  deformation,  another  equation,  describing  the  relationship  of 
shearing  forces  to  shearing  deformations,  is  required.  This  relation  is 

KA<l(|i  -  .  V  t  P  §  (51 


-O' 


wlujru  KAG  is  the  shear  rij'itlity  of  the  beam,  and  Ou/l)x  -  ip  is  tlic 
sliearing  slope..  Incorporating  Equations  4  and  5  into  tliu  equilibrium 
Equations  I  and  2  yields  the  following  Timoshenko  beam  equations  (see 
References  5,  11,  and  12): 


A(k-  i.)(K.r..+  m. 


N 

a  y  A(x  - 
"  1 


f.)  m.^. 

J  J  J 


(6) 


A  simplified  form  of  these  equations  can  be  found  in  many  stand¬ 
ard  texts  on  the  vibrations  of  continuous  systems  (see,  for  example. 
Reference  13).  The  provisions  for  attached  springs  and  masses  and 
the  inclusion  of  the  axial  compression  force  P  =  P(x)  are  not  usually 
given  explicitly  in  these  texts.  Inclusion  of  the  axial  compression  force 
in  the  equations  is  sometimes  important  in  missile  applications  because 
the  axial  acceleration  of  the  missile  in  flight  can  affect  the  bending 
dynamics. 

At  this  point  it  is  useful  to  review  briefly  what  is  meant  by  nor¬ 
mal  modes  and  describe  some  of  their  properties.  If  an  external  force 
acting  on  the  missile  varies  sinusoidally  at  a  single  frequency,  and  this 
frequency  is  slowly  varied,  we  would  observe  a  series  of  resonances  at 
a  number  of  discrete  frequencies.  At  these  frequencies,  the  amplitude 
of  vibration  becomes  quite  large,  and  each  point  of  the  missile  reaches 
its  maximum  displacement  at  the  same  time,  describing  a  definite  pat¬ 
tern  in  space.  This  pattern  or  mode  shape  is  commonly  called  an  eigen¬ 
function  or  normal  mode  (given  by  <t>j^(x),  here),  and  the  frequency  . 

at  which  it  occurs  is  called  an  eigenvalue.  A  very  important  property 
of  these  modes  is  that  the  motion  of  the  missile  in  any  one  mode  is  inde¬ 
pendent  of  the  motion  taking  place  in  any  of  the  others,  except  through 
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coupling  from  oxtiM'ital  (orctr.H  (sec  latc-r),  TMh  Ih  duo  to  tlio  ortlio- 
gouality  property  of  the  normal  modc*«,  and  bocaufli?  of  tliia  indcpondonco 
it  is  usually  OKpc.'dicnt  to  describe  the  missile  motion  in  terms  of  thosf! 
m  od  (!  a. 

Tlie  mathematical  expresHions  which  describe  tlic  ortliogonality 
of  the  normal  modes  can  be  obtained  directly  from  the  differential 
equations  of  motion.  The  derivation  of  these  orthogonality  conditions  is 
straightforward  and  follows  the  standard  procedure  given  in  many  texts. 
The  orthogonality  conditions  for  Equation  6  are;  (see  Reference  11) 

f  N  (M  ,m  =  n 

(pA<|.^<|>„  f  dx  t  Y  m  (4>  " 

'I  J-* 


I  ^  -  *„) 

r  d  d(j)  N 

-  I  r  K.^. 

"i  j=i  J  J" 


N  /M  cj  ,  m  =  n 

+  a  y  m.(i)j.  1.  +  U-.  t.  )  =(  ^ 

jWjn  VSna  \  ^  ^ 


where 

is  the  eigenfunction  for  deflection,  and 
i(i^  is  the  eigenfunction  for  bending  slope 

and  0  .  =0  (x)|  ,  ,  that  is,  0  (x)  evaluated  at  x  i  .,  and  f  .  is  the 

■^jm  .  m'  '(x=ij'  j'  '’jm 

displacement  of  the  jth  mass  in  the  mth  mode,  and  the  subscripts  m  and 

n  refer  to  two  different  modes.  The  constant  is  the  generalized  mass 

of  the  nth  mode,  cj  is  its  circular  frequency,  and  a  is  the  axial  accel- 
n  X 

eration  of  the  missile.  If  we  denote  the  amplitude  of  the  nth  mode  at  any 
instant  of  time  by  q^(t),  these  amplitudes  can  be  used  as  generalized 
coordinates. 


H(!ttinK  ii[i  tlir:  kiiiotir  Jliui  potoiitial  in  torma  of  tIuiBii  nn(0i 

UIkI  iiHinn  ortlio^jonality  conillllf)nH  7,  tlio  criuntioiis  of  iriotion  for  tlio  ^1,^(0 
can  bo  roduccd  to  the  following  form:  (huo  Reference  12  for  more  detail) 


M  fi  (t)+  M  2h  cj  q  t  M  u 

II  'll  ti  II  n  n  a  ii 


'q„(t)  .  Q,,(t) 


z{pA^) 


-  mjlaj  f> 

4  Q^(t)  +  5^(t)  (8) 


where  Qjj(t)  is  the  generalized  force  acting  on  the  nth  mode.  This  force 
can  be  found  in  the  usual  way  by  giving  the  nth  mode  a  virtual  displace¬ 
ment  Aq^  and  finding  the  work  done  by  all  the  external  forces  acting  on 
the  missile.  The  presence  of  only  one  mode  amplitude  q  in  Equation  8 
demonstrates  the  usefulness  of  the  bending  modes  being  independent  of 
each  other.  and  are  defined  later. 

To  return  to  the  problem  at  hand,  Equation{6)  can  be  solved,  sub¬ 
ject  to  appropriate  boundary  conditions,  to  yield  a  set  of  bending  modes 

Jjc 

(b  (x),  Ji  (x).  For  the  missile  in-flight,  the  ends  of  the  missile  are  free, 
'n  n 

and  the  engine  can  be  treated  as  an  integral  part  of  the  beam,  clamped 
with  zero  control  angle  6;  disconnected  from  the  missile  altogether;  or 
attached  to  the  missile  by  some  intermediate  constraint.  The  resulting 
equations,  showing  the  interaction  between  the  bending  and  the  control 
system,  differ  slightly  depending  on  what  representation  of  engine  attach¬ 
ment  is  made  for  the  computation  of  the  bending  modes.  The  equations 
that  follow  assume  the  modes  to  have  been  computed  with  the  engine 
clamped  at  6  =  0. 

* 

Some  practical  procedures  for  computing  the  bending  modes  are 
discussed  in  References  8  and  16. 
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Typical  frcc-frec  tniaailc  bending  modes  and  frequencies  for 
various  flight,  times  are  shown  in  Figure  3,  It  ia  noted  that  the  first 
free-free  mode  frequency  varies  by  less  than  10  percent  in  the  flight 
time  interval  from  launcl.  to  100  secends 

The  influence  of  shear  deformations,  rotatory  inertia,  and  axial 
acceleration  on  the  missile  bending  frequencies  is  shown  in  P’igure  4. 

In  the  first  mode  the  inclusion  of  shear  deflection,  noted  as  KAG  =  n 
and  the  exclusion  of  rotatory  inertia  noted  as  pi  =  0  increases  the  fre¬ 
quency  by  0.5  percent  over  the  frequency  if  both  effects  are  mciuded 
Note  that  the  ex'.:lusion  of  stiear  deiorrnatton  changes  the  first  mode 
frequency  by  12.  i  percent.  If  axial  acceleration  is  omitted  in  the  bend¬ 
ing  mode  analyses,  the  first  mode  frequency  is  increased  by  27  percent 
and  one  percent  for  constrained  tprelaunch)  and  free  flight  conditions 
respectively. 

2 . 2  Propellant  Dynamics  —  Sloshing 

The  description  of  the  missile  dynamics  for  missiles  with  liquid 
propellants  is  complicated  by  the  interaction  between  the  airframe  and 
the  propellants.  As  the  missile  bends  and  rotates,  the  liquid  propel¬ 
lants  move  back  and  forth  in  their  tanks  as  shown  in  Figure  5,  creating 
very  appreciable  forces  on  the  missile. 

It  can  be  shown  that  only  the  fundamental  (lowest  frequency)  mode 
of  propellant  motion  contributes  significantly  to  these  forces.  The  funda 
mental  mode  of  motion,  depicted  in  Figure  5  is  called  the  "sloshing" 
mode.  The  resonant  frequencies  of  these  sloshing  modes  are  usually 
quite  low  (approximately  0.5  cps)  and  they  interact  severely  with  the 
autopilot.  Furthermore,  since  the  sloshing  frequencies  are  consider¬ 
ably  below  the  first  free-free  "bending’  frequency  of  the  missile  struc¬ 
ture,  the  external  airframe  rotates  back  and  forth,  essentially  as  a 
rigid  body  in  some  of  the  "sloshing"  modes,  in  a  direction  opposite  from 
the  fluid  sloshing  in  order  to  preserve  angular  momentum.  This  means 
that  the  interaction  between  sloshing  and  aerodynamics  is  also  rather 
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Figure  3.  Typical  Free-Free  Missile  Bending  Modes 


Figure  4.  Influence  of  Shear  Deflection,  Rotatory  Inertia,  and  Axial  Acceleration 
on  Bending  Frequencies 


LIQUID  FREE  SURFACE 


Figure  5.  Fluid  Sloshing 


strong  since  the  entire  missile  picks  up  an  angle  of  attack  during  the 
sloshing  motion.  The  interaction  between  aerodynamic  foi'ces  and  those 
missile  modes  that  consist  mainly  of  structural  bending  is  not  nearly  as 
strong;  only  a  local  angle  of  attack  is  generated,  and  the  resulting  aero¬ 
dynamic  forces  are  usually  negligible. 

One  possible  scheme  for  eliminating  the  sloshing  problem  would  bo 
to  divide  the  propellant  tanks  into  small  compartments,  thereby  increas¬ 
ing  the  slosh  frequencies  and  reducing  the  sloshing  forces.  Unfortunately, 
because  of  weight  considerations,  such  a  procedure  is  generally  highly 
impractical.  However,  small  ring  baffles  are  usually  placed  at  various 
levels  in  the  tanks.  These  baffles  do  not  change  the  sloshing  frequencies 
but  do  provide  natural  damping  by  creating  turbulence  in  the  regions  of 
high  liquid  velocity. 


The  sopai'iile  airframe  and  propellant  dynamics  may  be  combined 
to  yield  what  are  known  as  bci^ding- sloshing  modes.  To  obtain  these 
combined  modes,  the  fluid  is  treated  as  concentrated  along  the  missile 
center  lino,  and  the  sloshing  modes  in  eacii  tank  are  represented  by 
attached  springs  and  masses.  The  spring-mass  systems  are  selected 
to  give  the  correct  slosh  frequencies  and  horizontal  force,  and  their  at¬ 
tachment  points  are  selected  to  give  the  same  moments  about  the  tank 
center  of  gravity  as  is  produced  by  the  fluid.  The  rotatory  inertia 
assigned  to  the  fluid,  pi  in  Equation  6,  is  adjusted  so  that  the  overall 
inertia  of  the  fluid  in  each  tank  is  the  same  as  that  found  by  solving  the 
fluid  dynamics  problem  of  a  tank  capped  at  the  quiescent  liquid  surface. 
Equation  6,  when  coupled  with  appropriate  equations  defining  K.,  will 
yield  combined  bending- sloshing  natural  frequencies  and  mode  shapes. 
Those  readers  who  wish  to  pursue  the  propellant  dynamics  problem 
further  will  find  a  discussion  of  the  sloshing  problem  in  References  7 
and  14. 

2.  3  Control  System 

The  control  system  or  autopilot  is  a  very  important  part  of  a  mis¬ 
sile  system.  The  thrusting  missile  is  unstable  in  flight  and  would  soon 
tumble  hopelessly  without  a  properly  designed  autopilot.  The  presence 
of  structural  bending  vibrations  makes  this  stability  problem  much  more 
complicated.  We  can  easily  demonstrate  this  for  ourselves  by  first 
taking  a  relatively  rigid  yardstick  in  the  hand  and  balancing  it  upright  by 
applying  only  vertical  and  lateral  forces  at  its  base.  The  yardstick  is 
balanced  with  ease.  If  this  yardstick  is  now  replaced  by  a  long  flexible 
reed  or  wire,  the  task  of  balancing  becomes  very  difficult. 

The  essential  function  of  the  autopilot  is  to  sense  the  missile 
attitude  and  angular  velocity  and  provide  engine  angular  deflections,  and 
hence  a  lateral  component  of  the  thrust,  in  accordance  with  these  sensed 
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quantities  so  as  to  maintain  the  missile  at  a  desired  attitude,  The 
problems  of  control  system  design  are  treated  extensively  elsewhere, 

BO  that  in  this  report  the  autopilot  will  appear  to  a  largo  extent  as  an 
unobtrusive  black  box.  In  order  to  provide  tlio  reader  with  something 
more  concrete,  however,  a  simplified  hydraulic  servo  will  be  presented 
in  block  diagram  form. 

2,  4  Propulsion  System  * 

The  propulsion  system  characteristics  will  enter  the  equations  of 
this  section  only  in  a  fairly  superficial  way.  The  propulsion  system  is 
taken  to  consist  of  nozzles  having  mass  and  moment  of  inertia  which 
rotate  about  gimbal  points  and  thereby  change  the  direction  of  a  steady 
thrust  vector. 


( 


( 
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3.  SIMPLIFIED  EQUATIONS  OF  MOTION 


The  equations  of  motion  that  follow  (as  in  reference  5)  assume 
that  all  slowly  varyinj.;  quantities  such  as  missile  mass,  axial  accel¬ 
eration,  liquid  levels,  mode  shapes  and  frequencies,  etc,,  can  he 
treated  as  fixed.  It  must  therefore  he  empiiasized  that  the  resulting 
equations  are  valid  only  for  a  few  seconds  at  any  selected  point  along 
the  trajectory;  however,  this  is  sufficient  to  study  '’high''  frequency 
in-flight  dynamics, 

3.  1  Missile  Deflections 

In  order  to  describe  the  motion  of  the  missile  in-flight,  it  is 
necessary  to  know  the  motion  cf  each  section  of  the  missile  as  a 
function  of  time.  Once  the  engine  start  transient  has  damped  out 
(this  problem  is  considered  in  the  next  section),  variations  in  axial 
forces  at  frequencies  high  enough  to  excite  longitudinal  vibrations 
(with  the  exception  of  engine  noise,  which  is  an  entirely  different 
type  of  problem)  are  very  small  so  that  all  longitudinal  deflections 
are  neglected  here.  The  position  of  each  section  of  the  missile  can 
be  described  by  two  functions,  u(x,  t),  the  lateral  deflection  of  each 
point  along  the  missile;  and  ip(x,  t),  the  angle  through  which  each 
section  rotates  relative  to  the  moving  axes.  See  Figure  6. 

These  functions  can  be  written  as 


oo 

u(x,  t)  =  C^(t)  fxC^(t)  +  Y. 

n=  1 

oo 

i|i(x,t)  -  Cj(t)  t  Y  ^n<^> 

n=  1 

where  x  is  measured  from  the  missile  center  of  gravity. 
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Figure  6.  Geometry  and  Coordinate  System 


The  fiinction8  0^(x)  and  are  the  aloahing-bending  modes 

described  earlier,  and  are  the  generalized  coordinates  of  these 

modes  and  represent  the  amplitude  that  must  be  assigned  to  each 
mode  at  any  instant.  Tlie  lime  functions  C^(t)  and  Cj(t)  represent  a 
rigid-body  translation  and  rotation,  respectively,  of  the  missile. 

These  quantities  can  be  defined  so  that  the  total  angular  and 
lateral  momentum  with  respect  to  the  moving  axes  is  independent  of 
the  engine  deflection  6.  To  do  this,  and  are  adjusted  so  that 
any  momentum  created  by  engine  motion  is  cancelled  by  rigid-body 
motions  of  the  missile  as  a  whole.  Applying  this  condition,  and 
C  j  are  given  by 

tTO  Z  A 

C  - - It—  6  =  -K  6 

O  M  O 

=  1  =  r- (’e  +  <>'> 

Where  is  the  mass  moment  of  inertia  of  the  entire  missile  (including 
engines)  about  the  missile  center-of-gravity,  is  the  mass  moment 
of  inertia  of  the  swiveling  engines  about  their  gimbal  points,  and  the 
other  parameters  are  shown  in  Figure  6. 

With  these  definitions  of  and  Cj,  the  total  angular  momentum 
depends  only  on  and  0.  The  contribution  of  is  found  by  first 
writing  the  equation  of  moment  equilibrium  between  the  inertia  forces 
and  moments  and  the  moment  due  to  the  engine  force  for  the  missile 
deflected  into  an  extreme  position  in  the  nth  mode: 


(12) 
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where  T  and  T  are  enuine  thrusts,  described  in  more  detail  later. 
After  multiplying  through  this  equation  by  inspection  of  the  first 
two  terms  shows  them  to  be  the  moment  of  momentum  about  the  missile 
center  of  gravity  due  to  a  coordinate  velocity  q^^.  This  equation  can 
then  be  used  to  determine  the  following  simple  expression  for  this 
momentum  in  the  positive  0  direction: 

H  =  [t  +  T  ]  --  q  =  K  q  (13) 

n  \  a  cj  2  2n  m 

'  u 

n 


The  contribution  to  the  angular  momentum  due  to  0  is  merely 


3.  2  Rigid-Body  Equations 

The  external  forces  that  must  be  considered  in  writing  the  rigid- 
body  equations  are  shown  in  Figure  7. 

Since  there  may  be  more  than  one  engine,  and  not  all  the  engines 

may  swivel  in  each  plane,  the  total  thrust  is  broken  into  two  components; 

T  is  the  total  thrust  of  the  swiveling  engines,  and  T  is  the  total 
c  s 

thrust  of  the  engines  stationary  with  respect  to  the  tangent  to  the  missile 
axis  at  the  attach  point. 


The  equation  governing  rotation  (pitching)  about  tlie  missile  center 
of  gravity  is  obtained  by  taking  moments  about  the  center  of  gravity  with 
the  result  (see  reference  12) 

CD 


y'  K,  q  +  I  9  -  ~a  f  xw(x)  dx  +  f  T  (6  -  u  -  T  u  '  1  I 
2nn  c  T  Lev  ej  seJe 

-(t  +  T  ]u 
V  c  3/  e 


(14) 
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Figvire  7.  Forces  Acting  on  Missile 
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Here 


=  p  -  0  +  E  ^  +  Kj  6  -  ^  (15) 

n=  1  o  o 

is  the  total  angle  of  attack  producing  aerodynamic  forces.  The  velocity 
is  the  component  of  atmospheric  wind  velocity  perpendicular  to  tlie 
x-axis.  The  angle  of  attack  due  to  the  sloshing  is  represented  by 
a^q^(t)  where  is  the  angle  through  which  the  missile  airframe  rotates 
per  unit  of  nth  "slosh  mode"  (a  bending  -  sloshing  mode  which  consists 
mainly  of  sloshing)  amplitude.  The  summation  in  Equation  (14)  is 
taken  only  over  N  such  "slosh  modes,"  It  is  implicitly  assumed  in 
Equation  (14)  that  the  aerodynamic  force  depends  linearly  on  the  angle 
of  attack,  and  that  w(x)  is  the  running  load  per  unit  angle  of  attack.  The 
dependence  of  the  aerodynamic  forces  on  higher  powers  of  a,  and  on 
other  quantities,  such  as  the  angular  velocity  of  the  missile,  are  neglected 
here  for  simplicity.  These  additional  forces  are  usually  small,  and  for 
many  studies  it  is  legitimate  to  neglect  them. 

The  angle  u^'  (the  prime  indicates  differentiation  with  respect  to 
x)  is  the  slope  of  the  missile  axis  at  the  engine  gimbal  point.  It  is 
obtained  by  differentiating  Equation  (9)  and  evaluating  at  x  = 

oo 

u  '(t)  =  K  6  t  y  9  '  (-f„iq  (In) 

The  equation  governing  rigid-body  translation  is  readily  obtained  by 
summing  forces  in  the  y  direction  (Figure  7)  with  the  result 

My  =  ■a.p  w(x)  dx  +  Mg  0  +  (  +  T^)u^'  -  T^6  (17) 

where  M  is  the  total  mass  of  the  missile  and  g  is  the  acceleration 
caused  by  gravity. 


3,  3  Miaflile  Bending  EquatlonB 

The  bending  equations  (Including  sloshing)  are  most  easily  derived 

by  first  giving  each  mode  a  virtual  displacement  and  observing 

the  work  done  by  the  forces  acting  on  the  missile.  If  this  work  is 

divided  by  the  virtual  displacement  we  obtain  the  generalized  force 

Q  to  be  used  in 
n 


-t-  2b  u  q 
n  n^n 


+  CO  q 


n 


Q 

n 


-  w~ 

n 


(18) 


where  is  the  nth  generalized  mass  defined  by  Equation  7  and  (5^  is 
defined  in  Equation  8.  The  generalized  force,  Q^,  is 


Q 


n 


w(x)  ^jj{x)  dx 


Ft  [s  -  u  ') 

-  T  u  ^  (-1  ] 

L  C  1  e) 

s  e  J  ’^n  \  ej 

(19) 


In  Equation  (18),u^  is  the  circular  frequency  of  the  nth  mode  and  is 
the  damping  ratio.  Damping  ratio  b^  is  either  determined  experiment¬ 
ally  by  shaking  an  assembled  missile,  or  assumed  conservatively  small 
(0.  5  to  1  percent  of  critical).  It  should  be  mentioned  that  the  damping 
of  a  structure  can  seldom  be  represented  in  this  way.  In  order  that 
this  representation  be  strictly  correct,  all  damping  forces  would  have 
to  be  of  the  viscous  type  (that  is,  proportional  to  velocity)  and  would 
have  to  be  distributed  along  the  length  of  the  missile  in  the  same  way  in 
which  the  mass  and  rotatory  inertia  are  distributed,  or  in  the  same  way 
in  which  an  appropriate  combination  of  the  elastic  forces  are  distributed. 
If  the  damping  forces  are  distributed  in  any  other  way,  application  of 
Lagrange's  eqviations  would  result  in  the  equations  of  motion  for  the 
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bcndlnf(  modoH  having  velocity  coupling  tltrough  the  viacoua  damping 

forccH,  That  la,  each  equation  correaponding  to  Equation  (18)  would 

have  additional  terma  of  the  form  C  q  ,  where  m  =  1,2 .  For 

m  m 

reasonably  small  values  of  damping,  however,  (corresponding  to  a 
few  percent  of  critical  damping  in  any  mode),  these  coupling  terma 
have  a  negligibly  small  effect  on  the  motion  (reference  15)  and  can  be 
omitted  as  waa  done  in  Equation  (18).  Similarly,  if  the  damping  is  of 
the  hysteresis  type  rather  than  viscous,  it  can  still  be  represented  by 
an  equivalent  viscous  damping  ratio  b^  for  sufficiently  small  damping. 

3.  4  Control  System  and  Engine-Swiveling  Equations 

The  control  system  will  be  assumed  to  sense  two  quantities  from 
gyros  mounted  on  the  missile  airframe.  These  quantities  are  the 
attitude  of  the  missile  at  the  gyro  station,  given  by  =  9  +  u  ',  and 

the  angular  velocity  at  that  station,  given  by  0— where 

f  .D 

CX> 

=  (201 


If  the  missile  begins  to  go  off  course  by  building  up  an  attitude  0,  or 
if  the  gyros  think  the  missile  is  going  off  course  because  the  missile 
bends,  giving  rise  to  a  u^',  the  control  system  calls  for  a  corrective 
engine  angle,  that  is,  a  lateral  thrust  component.  The  engine  command 
angle  6^  is  formed  as  a  linear  combination  of  the  difference  between 
a  desired  attitude  9^  and  the  sensed  attitude  9p,g  and  the  sensed  rate 
0pg.  This  is  expressed  by 

= '^oK-^fb)  - Vfb  <^'1 
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where  Kg  and  are  gain  constants.  This  engine  command  signal  Is 
modified  by  compensation  n:;tworks  which  shape  the  phase  and  amplitude 
of  the  signal  to  give  satisfactory  stability.  The  resulting  signal  6  is 
used  to  open  an  hydraulic  valve,  as  shown  in  Figure  8. 


Figure  8.  Schematic  of  Autopilot  Actuator 

The  fluid  controlled  by  this  valve  flows  into  an  hydraulic  actuator 
which  is  connected  directly  to  the  engine  bell.  The  position  of  the  actuator 
rod  6  is  sensed  and  subtracted  from  the  compensated  engine  command 
angle  6^  so  that  it  is  actually  the  error  signal  5^  -  5^  which  determines 
the  rate  of  flow  in  the  valve.  If  the  engine  angle,  as  measured  by  the  posi¬ 
tion  of  the  actuator  rod,  is  the  same  as  the  command  angle  6^,  no  further 
fluid  flows  and  the  actuator  does  not  move. 
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Going  one  step  farther,  if  the  compressibility  of  the  hydraulic 
fluid  is  taken  into  account,  there  is  an  apparent  compressibility  flow 
which  must  be  subtracted  from  the  flow  initiated  by  the  error  signal 
6,.  -  6,.  This  entire  process  is  represented  conveniently  in  block 
diagram  form  in  Figure  9.  The  variable  s,  used  in  this  diagram,  is 
the  usual  Laplace  operator. 

From  Figure  9  it  is  clear  that  the  control  flow  in  the  valve  is  pro¬ 
portional  to  except  for  a  small  delay  or  valve  action  time  t.  From 
this  flow  rate  the  compressibility  flow  is  subtracted  and  the  result¬ 
ant  flow  is  divided  by  Ar  (defined  by  Figure  8)  and  integrated  (multiplied 
by  l/s)  to  give  the  actuator  position  6  ,  as  shown.  The  difference 
between  the  engine  position  6  and  the  actuator  position  6  is  multiplied 
by  the  equivalent  torsional  spring  constant  of  the  actuator  to  give 
the  torque  acting  on  the  engine  bell.  Following  around  the  hydraulic 
loop,  this  torque  divided  by  Ar  gives  the  pressure  in  the  hydraulic 

Cylinder.  This  pressure,  multiplied  by  the  effective  volume  V  under 

o 

compression  and  divided  by  the  bulk  modulus  of  elasticity  b  of  the 
fluid,  gives  the  change  in  volume  of  the  fluid  caused  by  the  pressure. 

The  rate  of  change  (multiply  by  s)  of  volume  is  the  compressibility  flow 
Q^.  The  torque  acting  on  the  engine  bell  causes  it  to  move  according  to 

I.‘  +  B,S  +  K^5  =  K^(6^  -  s).  (22) 

where  is  the  torsional  spring  constant  between  the  engine  bell  and 
the  missile  airframe  (caused  by  bellows,  pendulum  action  resulting 
from  axial  acceleration,  etc,  )  and  represents  any  viscous  damping 
present.  Equation  (22)  is  given  in  block  diagram  form  in  Figure  9. 

This  equation  is  the  only  equation  in  this  development  of  the  integrated 
dynamic  equations  which  is  basically  incomplete.  Many  terms  have 
been  omitted  from  the  other  equations,  but  they  were  omitted  because 
they  were  small.  In  the  case  of  Equation  (22),  however,  all  the  inter¬ 
actions  between  engine  swiveling,  bending,  and  rigid-body  dyanmics 
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Figure  9.  Block  Diagram  of  Autopilot  and  Engine  Dynamics 


have  been  neglected  for  the  sake  of  expediency.  These  interactions  are 
important  in  relation  to  other  effects  included  in  the  analysis,  but  most 
of  them  arise  in  a  manner  so  subtle  that  Lagrange's  equations  must  be 
resorted  to  to  demonstrate  their  existence.  Such  a  procedure  is 
lengthy  but  straightforward,  and  the  enterprising  reader  can  derive  the 
more  complete  engine-swiveling  equation  given  by  (see  reference  12) 


-Kmz-K,(l  +miz]l6  +  mza  (6-u'^ 
oe  lie  e  c  JJ  e  x\  e  J 


,  - .  .  _  “  . .  f-b 

\ 

I  +mfzl0tmza  t  )  q  /  p 

\  e  e  e  J  e  y  Jo  ^ 

'AV  -  PHJ 

=  K,(5^  -  6)  -  bJ  -  Kj6  (Z31 

where  b  is  the  distance  from  the  engine  gimbal  to  the  nozzle  exit  plane, 

a  and  a  are  the  rigid-body  accelerations  in  the  x  and  y  directions, 

X  y 

and  is  the  spring  constant  between  the  engine  bell  caused  by  a 
direct  spring  action,  such  as  from  a  bellows.  In  spite  of  its  being  in¬ 
complete,  the  use  of  Equation  (22)  in  conjunction  with  the  other  equa¬ 
tions  developed  here  would  generally  give  adequate  engineering  inform¬ 
ation  concerning  in-flight  dynamics. 

3.  5  Block  Diagrams 

Following  the  procedure  used  for  the  control  equations,  we  can 
convert  the  bending  and  rigid-body  equations  into  Laplace  notation  and 
draw  a  block  diagram  which  represents  these  equations  completely. 

Such  a  diagram  is  given  in  Figure  10.  Only  one  block  representing  a 
typical  bending-sloshing  mode  is  shown  in  this  diagram;  however,  it  is 
understood  that  there  are  as  many  such  blocks  as  there  are  modes 
being  used.  In  most  cases  inclusion  of  the  "slosh  modes"  (one  for  each 
tank)  and  two  or  three  bending  modes  is  sufficient.  The  summation 
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Figure  10.  Block  Diagram  Showing  Detailed  In-Flight  Interactions 


g;vfcn  on  tho  patJis  loading  out  of  thu  bending- wloahing  block 
indicate  that  a  auinination  Is  to  be  made  over  the  modes  of  the  quanti¬ 
ties  in  question. 

A  itudy  of  Figure*  10  is  enlightening  in  that  tlu;  flow  of  cause  and 
effect  of  each  subsystem  on  the  other  can  be  seen  without  having  to 
interpret  several  equations  simultaneously.  In  addition,  several 
significant  quantities  such  as  the  engine  lateral  force  caused  by  bend- 

I 

ing  and  the  total  angle  of  attack  q.^,  I  see  Equation  ( 1 5)  j  and  their  effect 
on  the  system  become  apparent.  The  complexity  of  Figure  10  detracts 
somewhat  from  its  usefulness  if  we  do  not  care  to  follow  all  the  details; 
therefore,  a  functional  diagram  showing  only  the  directions  of  the  inter¬ 
actions  is  given  in  Figure  11. 

3.  6  Use  of  Equations 

The  in-flight  dynamics  equations  developed  herein  have  many  uses. 
From  the  standpoint  of  the  control  system  engineer,  their  most  impor¬ 
tant  use  is  to  allow  him  to  study  the  stability  of  the  missile  in-flight 
and  to  indicate  the  form  of  autopilot  compensation  required  to  obtain 
stability.  The  structural  dynamicist  sees  them  as  a  tool  to  determine 
the  dynamic  loads  on  critical  sections  of  the  missile.  These  loads 
arise  from  autopilot  limit  cycles,  discrete  and  random  gusts,  staging 
transients,  and  other  rapid  maneuvers  where  the  action  time  is  of  the 
order  of  the  period  of  one  of  the  bending  modes. 

For  the  control  engineer  these  equations  are  complete  as  they  are 
(with  the  exception  of  the  simplifications),  but  the  structural  dyna¬ 
micist  must  modify  them  slightly  to  insure  satisfactory  convergence  of 
modal  solutions,  that  is,  a  series  expansion  in  terms  of  normal  modes. 
This  happens  because  the  control  engineer  is  interested  only  in  dis¬ 
placements,  whereas  the  structural  dynamicist  is  interested  in  bending 
moments  and  shear  forces  as  well.  The  convergence  of  series  expan¬ 
sions  of  bending  moments  and  shear  forces  is  very  poor  for  a  beam 
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Figure  11.  Functional  Block  Diagram 


( 


with  a  concunt  rated  load  (the  engino  force),  and  even  in  the  ab.<ienco  of 
any  dynamic  effects  at  all,  many  inodes  would  be  needed  to  obtain 
satisfactory  accuracy.  On  the  other  hand,  series  expansions  of  dis¬ 
placement  and  slope  converge  quite  rapidly. 

The  difficulties  associated  with  computing  bending  moments  and 
shear  force  can  be  minimized  by  employing  the  so-called  "mode 
acceleration"  technique  (reference  16).  The  essential  feature  of  this 
method  is  that  it  separates  the  pseudostatic  loads  (loads  found  by  rigid- 
body  techniques)  from  the  loads  arising  from  the  dynamic  response  of 
the  missile.  Thus,  the  poor  convergence  caused  by  discontinuities  in 
the  loading,  such  as  a  concentrated  force,  are  taken  into  the  pseudostatic 
load,  and  the  additional  load  caused  by  dynamic  bending  usually  con¬ 
verges  much  faster.  For  example,  the  bending  moment  at  any  station, 
along  a  missile,  can  be  written  as 


oo  q  (t) 

M(x,  t)  =  M  (x,  t)  -  Y. 

s  ,  n  c 

n=  1  w 

n 


where  M^(x,  t)  is  the  pseudostatic  load  neglecting  dynamic  bending,  and 
m^(x)  is  the  nth  eigenfxmction  representing  the  bending  moment  at  any 
station.  The  same  moment  expressed  without  the  use  of  the  mode 
acceleration  procedure  is  given  by 


oo 

M(x.t)  =  y  m  (x)q  (t)  (Z5) 

,  n  n 
n=  1 

which  converges  more  slowly  because  of  the  absence  of  in  the 

denominator.  Equation  (24)  assumes  small  damping  and  can  easily  be 
derived  from  Equation  (25)  by  using 


q  +  u  q 
’n  n  ^n 


M 


(26) 


( 
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Solvjnjj  for  from  Equation  (26)  and  aubstitutinj^  into  Equation  (25), 
we  obtain 

GO  m  (x)  Q  (t) 

M(x,t)  =  Z  ^ 

n-  1  Co  M 
n  n 

But 

CO  m  (x)  Q  (t) 

Zn  n 

2 

n=  1  M  (o 
n  n 


^  m^(x)  — ^ 

n=  1  CO 


is  simply  the  moment  that  results  if  the  system  of  forces,  including 
rigid-body  inertia  forces,  giving  rise  to  Qj^(t)  is  applied  "statically" 
since  it  comes  from  the  solution  of  Equation  (26)  for  q^  =  0.  Thus 
Equation  (27)  is  identical  with  Equation  (24).  A  term 


oo  2b  q 

y  m  (x)  — ^ 
^ ,  n'  CO 
n=l  n 


would  appear  on  the  right-hand  side  of  Equation  (24)  if  damping  were 
present  in  Equation  (26),  but  the  damping  ratios  b^  are  usually  quite 
small,  and  these  terms  can  be  neglected. 
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4.  .ILLUSTRATIVE  EXAMPLES 


4.  1  Limit  Cycle  Response 

The  interaction  between  control  ayatem  dynamics  and  air/rame 
bending  dynamics  can  force  the  control  system  into  a  nonlinear  mode 
of  operation  which  results  in  a  limit  cycle.  In  such  a  situation  the 
control  system  causes  the  engines  to  gimbal  continuously,  searching 
for  the  correct  trim  condition,  that  is,  the  control  system  "hunts"  for 
the  correct  engine  swivel  angle  6.  The  data  in  Figure  12  illustrate 
the  effects  that  a  control  system  limit  cycle  can  have  on  dynamic  loads 
The  data  were  taken  from  an  actual  case  (time  of  flight  corresponds  to 
maximum  dynamic  pressure)  in  which  the  damping  in  the  control  sys¬ 
tem  hydraulics  became  negative  for  large  engine  angular  velocities  6. 
The  negative  damping  caused  a  limit  cycle  whose  frequency  was  near 
the  first  bending  frequency  of  the  missile.  The  control  system  had  to 
be  redesigned  to  eliminate  the  limit  cycle. 
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Figure  12.  Dynamic  Loads  Generated  by  a  Control  System  Limit  Cycle 


4.  2  Reaponae  Dm;  to  Gnat 


{ 


It  is  recognized  that  missile  inflight  response  analyses  are 
complex.  This  is  due  to  the  interactions  of  the  control  system,  aero¬ 
dynamics,  airframe  bending  dynamics,  and  the  engine  swiveling  dy¬ 
namics  as  illustrated  in  Figure  11.  However,  such  analyses  are 
tractable  on  electronic  computers.  Typical  results  for  the  inflight 
response  of  a  ballistic  missile  during  its  flight  through  a  triangular 
shaped  (A)  gust  are  shown  in  Figures  13  and  14,  This  analysis  was 
performed  on  an  analog  computer.  In  Figures  13  and  14  the  following 
nomenclature  is  adopted; 


"N 


generalized  bending  acceleration  in  the  nth  mode 
whe re  n  -  1,2,3. 

angle  of  attack  representing  the  gust,  used  as  an 
input  or  driving  function. 

total  angle  of  attack, 
nozzle  vector  angle. 

psuedostatic  bending  moment  at  station  x, 
dynamic  bending  moment  at  station  x. 
total  bending  moment  at  station  x. 


A  "mode  acceleration"  technique  (as  discussed  previously)  is 
used  in  the  analysis  and  consequently  one  is  interested  in  the  generalized 
bending  accelerations  (q^).  recalling  that. 


M.j,(x,  t)  =  Mj^(x,  t) 


q„(t) 


n  - 1 


Mplx,  t) 
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Figure  14.  Inflight  Response  Due  to  A  Gust 
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Fioiii  Figuif  1  i  It  Id  clear  lliat  llie  iiiaj‘>r  contribution  to  M  i.omea 

r 

from  the  first  mode.  In  Figure  14  Mp  is  less  th;in  5  percent  of 
at  any  time,  consequently  id  approximately  equal  to  This 

result  13  typical  for  such  problems  if  instabilities  are  not  present  even 
if  other  shapes  of  gust  velonties  are  considered  such  as  random  or 
”  1  -cos "  gusts. 


5.  PRELAQNCH-POSTLAUNCH  DYNAMIC  LOADS 


In  Section  2  the  equations  necessary  for  the  computation  of  in-flight 
dynamic  loads  were  developed.  The  important  elements  of  the  in-flight 
loads  problem  are  the  control  system,  the  airframe  bending  dynamics, 
and  the  in-flight  aerodynamics.  In  this  section,  three  types  of  prob¬ 
lems  are  considered.  The  first  is  concerned  with  the  interactions  be¬ 
tween  airframe  bending  dynamics  and  ground  wind  aerodynamic  forces . 
The  second  type  of  problem  has  to  do  with  the  interactions  between  the 
axial  or  (longitudinal)  airframe  dynamics  and  the  propulsion  system. 

The  third  type  of  problem  is  concerned  with  the  response  of  a  con¬ 
strained  missile  due  to  ground  shock. 

These  problems  fall  into  a  general  class  known  as  prelaunch, 
postlaunch. problems,  depending  on  whether  the  dynamic  loads  occur 
before  or  as  a  result  of  releasing  the  missile. 

5. 1  Wind-Induced  Oscillations 

5.1.1  Theoretical  Considerations.  A  missile  in  the  prelaunch  condi¬ 
tion,  that  is,  erected  on  a  launch  pad  (Figure  15),  may  be  exposed  to 
ground  winds  which  give  rise  to  flow  around  the  missile  of  large 
Reynolds  number  R  (R-  10^  based  on  the  missile  diameter).  In  such 
a  flow  the  boundary  layer  on  the  missile  is  turbulent  and,  as  a  con¬ 
sequence,  random  pressure  fluctuations  act  on  the  missile.  These 
random  pressure  fluctuations  create  random  lift  and  drag  forces  which 
in  turn  can  induce  large,  oscillating  displacements  and  bending  mo¬ 
ments  in  the  missile. 


The  Reynolds  number  R  is  defined  to  be  R  =  Ud/u  ,  where  U  is 
the  undisturbed  velocity  of  flow,  d  is  the  missile  diameter,  and 
u  is  the  kinematic  viscosity  of  the  flowing  liquid. 
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Figure  15.  Missile  Exposed  to  Ground  Winds 


In  I'tri  ;ilr  flow  porpondlcvila r  to  a  circular  cylinder,  such  as  a 
missile,  a  groat  variety  of  changes  occur  with  an  increasing  Reynolds 
number.  Details  of  these  changes  can  be  found  in  Reference  17  and 
other  texts  in  fluid  mechanics.  Briefly,  if  R  is  in  the  range  of  40  to 
150,  the  "shedding"  of  vortices  is  regular.  The  eddying  motion  in  the 
wake  is  periodic  both  in  space  and  time,  and  the  flow  can  be  approxi¬ 
mated  by  the  well-known  Karman  vortex  street.  The  range  of  R 
between  150  and  300  is  a  transition  range,  in  which  the  vortex  shedding 
is  no  longer  regular.  For  R=*300,  the  vortex  shedding  is  "irregular." 

A  predominant  frequency  can  be  determined,  but  the  amplitude  is  ran- 

5 

dom.  Finally,  at  R  of  order  10  ,  the  separation  point  of  the  boundary 
layer  moves  rearward  on  the  cylinder  and  the  drag  coefficient  de¬ 
creases  appreciably.  The  flow  in  the  wake  at  these  large  Reynolds 
numbers  becomes  so  turbulent  that  the  vortex  street  pattern  is  no 
longer  recognizable.  Experimental  results  for  flow  around  cylinders 

5 

when  R  is  larger  than  10  are  presented  in  Reference  18. 

Consider  the  lateral  motion  of  the  missile  in  Figure  15.  The 
applied  lift  force  per  unit  length  along  the  missile  can  be  written  as 

L(x.  t)  =  q(x)d(x)Cj^(t)  (28) 

where  q(x)  is  the  dynamic  pressure,  d(x)  is  the  diameter,  and 
Cj^(t)  is  the  lift  coefficient  which  varies  randomly  with  time.  Actually, 

Cj^  is  a  function  of  Reynolds  number  as  well  as  of  time.  Therefore 
C^,  at  each  missile  station,  depends  on  the  diameter  of  the  missile  at 
that  station.  The  dependence  of  on  R  and  t  is  such  that  it  is  not 

possible,  in  general,  to  express  as  a  product  or  sum  of  a  function 

of  R  and  a  function  of  t.  However,  the  experimental  data  of  Reference 
18  indicates  that  C.  is  not  a  strong  function  of  R,  at  least  for 

5 

limited  ranges  of  R  in  excess  of  5  x  10  .  Therefore,  it  is  reasonable 
to  assume  that  a  value  of  Cj^,  associated  with  a  Reynolds  number 
based  on  a  representative  missile  diameter,  should  be  adequate  for 
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Equation  28.  This  conclusion,  being  founded  on  the  data  of  Reference 
18,  is  only  valid  for  two-dimensional  flow.  The  three-dimensional 
flow  characteristics  around  the  tip  of  a  missile  have  a  substantial 
effect  on  the  total  lift  and  drag  forces.  (See  References  19  and  20.) 
However,  in  the  absence  of  detailed  wind  tunnel  input  forces  on 
particular  missiles,  the  procedure  developed  herein  is  the  best  avail¬ 
able.  On  the  basis  of  limited  experimental  data,  it  can  be  anticipated 
that  the  results  obtained  will  be  conservative. 

The  lift  force  of  Equation  28  is  a  random  function  of  time  since 
C^{t)  is  random  in  time.  Problems  involving  random  time  varia¬ 
tions  are  usually  more  conveniently  treated  in  the  frequency  domain 
than  in  the  time  domain.  The  frequency  representation  of  a  random 
function  is  known  as  its  power  spectral  density.  Figure  16  is  a  typical 
power  spectral  density  of  lift  force  on  a  uniform  cylinder.  The  fre¬ 
quency  parameter,  S  =  Gjd/2TrU,  used  in  Figure  16  is  known  as  the 
Strouhal  number:  the  curve  is  normalized  such  that 

r“  ,  f“^/'ud'\  ,  2ttU 

Jo  ®^<S)dS  =  1  or  F^j  du  - 

By  definition,  if  P{cj)  is  the  power  spectral  density  of  the 
lift  coefficient  Cj^(t),  the  mean  square  value  of  is  given  by 


P(w)  du) 


By  comparing  the  latter  two  equations,  the  normalizing  factor 
K  =  P(u)/F(S)  is  evident  and 


PJw) 


■HU 


F(s) 


(29) 


(30) 
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Figure  16.  Normalized  Power  Spectrum  for  the  Lift  Force 
and  Drag  Force  at  Reynolds  Number  1.39x10° 
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where 


^  F(x)^^(x)dx  and  F(x)  =  q(x)d(x) 


For  such  a  forcing  function  it  is  possible  to  find  an  ex¬ 
tremely  simple  expression  for  the  root  mean  square  (rms)  of  the 
time  response  of  the  missile  by  assuming  that  the  damping  ratios 
b^  are  small.  The  function  T(t)  is  assumed  to  be  random  and 
expressible  in  terms  of  a  power  spectral  density  P(u). 

To  obtain  the  response,  first  notice  that  if  T(t)  equals 
a  coso,  the  steady-state  solution  of  the  equation  of  motion  is 


%  = 


~T 


aW  cos  (uit  +  0  ) 
n  n 


M  <i> 


n  n  \ 


^'ii . 


'.L 


+  2b 


(35) 


The  corresponding  displacement  is 


^  W  0  (x)  cos  (ut  +  0  ) 

u(x.  t)  =  ^  ^  ^ 


-  rr  r  \2lv  ^ 

M  w  {  1  -  —  +  2b  — 

n  n  \  /  \  n  u) 

I  V  V 


?r7T 


The  mean  square  response  (which  will  be  needed  as  explained  later) 
is  defined  to  be 


(36) 


2,  .  u 


t+2Tr/u> 


u  (x,  t)  dr 


(37) 


u^(x) 


wa 

2ir 


•t'(-2w/<. 


V 


W  W  p  tp  cos  (wT+  0  )cos(wT+  0  ) 
n  n  d  7“ 


n  m 


(38) 
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This  equation  can  be  integrated  by  term  to  yield 


-T—  2  V  ’  W  W  ^4 

u  (x)  /  )  - 2 - 

n  m  n  m  n  m 


W  W  6  6  cos  (9  -  9  ) 

n  m'^n'^m  '  n  m 


,(cj,  n)  ,(lj,  m) 


Equation  39  is  the  mean  square  response  to  a  force  distribu¬ 
tion  fixed  in  space  and  varying  sinusoidally  with  time  at  frequency  u. 
If  T(t)  is  considered  as  a  concentrated  force  acting  on  a  unit  dashpot; 
the  average  power  dissipated  is  a^/2,  that  is 

E  =  fi  'TV)dT  =  V 

Similarly,  if  T(t)  is  a  randorn  force  acting  on  a  unit  dashpot,  a 
power  spectral  density  P(i»>)  can  be  defined  for  this  force  (with 
certain  restrictions)  such  that  the  average  power  dissipated  by  fre¬ 
quency  components  in  the  range  w,  u  +  dco  is  P(w)  du.  Thus,  if  a 
random  force  of  the  form  given  by  f(x,  t)  acts  on  the  missile,  its 
response  can  be  obtained  by  using  Equation  39  and  superimposing 
the  responses  from  each  frequency  component  of  the  random  force, 
that  is,  replace  a  /2  by  P(w)duj  and  integrate  over  all  frequencies; 


—5 —  (-CD-.--  ■  W  W  6  p  cos  (9  -  9  )  P(uj) 

2,  ,  (  '  n  m^n’^m  '  n  m'  '  , 

u  (x)  =  I  >  - 2 - 2"^- - r-p^y- - -JJ2-  du  . 

J  0  /  ^  M  M  tj  w  '(uj,  n)'  ,(cj,  m)> 

n  m  n  m  '  '  ' 


n  m  n  m 


The  manipulations  which  led  to  Equation  40  have  been  essentially 
intuitive  and  without  mathematical  justification.  However,  a  mathe¬ 
matically  rigorous  treatment  of  the  same  problem  yields  exactly 
Equation  40. 
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It  is  the  nature  of  random  processes  that  only  a  statistical 
description  is  meaningful.  It  is  quite  often  impossible  to  con¬ 
struct  the  statistics  needed  for  a  complete  description;  however, 
the  mean  square  value,  which  is  the  second  statistical  moment,  is 
usually  obtainable.  If  the  statistics  of  a  random  process  are 
gaussian,  the  mean  square  value  provides  a  complete  description. 

It  is  for  these  reasons  that  the  particular  manipulations  leading  to 
Equation  40  were  undertaken. 

In  evaluating  the  integral  in  Equation  40,  a  simplification 

can  be  made  if  the  damping  ratios  b^  are  small  and  if  P(<j)  is  a 

smooth  function  as  compared  to  the  resonant  peaks  of  the  structure. 

In  such  a  case,  the  largest  contributions  to  the  integrand  are  those 

terms  for  which  n  =  m  because  near  each  resonance  u  u  the 

n 

radicals  in  the  denominator  become  very  small.  By  neglecting  the 
terms  for  which  n  ^  m. 


~Z 

u 


(x) 


d(j 


or 


2  2 

\  '  W  (4  .-oo 


-  W  =>  -vv  I 


P(u) 


M  “w  ‘  '0 

n  n  n 


/  '2 

2 

+  I  2b  — 

1  ^  / 

\  n  u) 

V  "/ 

\  n  y 

du 


(41) 


If  P|w)  is  a  constant,  no  further  assumptions  need  be  made 

because  each  term  in  the  foregoing  equation  is  easily  integrated. 

Remembering  that  the  integrand  is  a  function  with  very  steep  peaks 

at  each  w  ,  we  can  make  another  simplifying  assumption.  Since 
n 

the  integrand  has  these  steep  peaks  near  w  -  the  main  contri¬ 
bution  to  the  integral  comes  when  P((j)  -  P(u>^):  thus  a  good  approxi¬ 
mation  results  if  P(w)  is  taken  as  -  constant  for  each  term. 
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Then 


— 5 —  ’  W  ^(x)  ttw 

Z,  ,  n  ’  n  „/  I 

u  (x)  ^  ;  -7^"  "^b" 

- -  M  ui  n 

n  nr 


A  similar  expression  is  derived  in  Referen.e  Zl.  Substituting  foi 
P(oj^)  from  Equation  30  yields 


u  (x)  = 


dC^^  s 
8U  , 


W  ^(p  ^(x) 


.  b  M  ^ 
n  n  n  n 


It  follows  from  Equation  4  that  the  mean  square  bending 
moment  along  the  missile  in  the  lift  plane  is 


- y  T  r  di|i 

- 3 -  dC.  --W  EI(x)-^ 

.  .Z,  V  L  n  L  dx 

^  =  8U  /  o  w  2  3 

L _  b  M  w 

n  n  n  n 


F(S„) 


The  mean,  square  response  in  the  drag  plane  can  be  obtained  from  this 
expression  if  C  is  replaced  by  • 

The  other  wind  force  acting  per  unit  length  along  the  missile 
is  the  steady- state  drag  force  defined  as 


=  C_-q(x)d(x) 


where  =  steady-state  drag  coefficient 

q(x)  =  dynamic  air  pressure  at  station  x 

For  purposes  of  analysis,  the  responses  in  the  drag  and 
lift  planes  are  added  vectorally.  For  example,  the  peak  displace¬ 
ment  of  the  missile  at  station  x  is  defined  as 
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1/2 


(40) 


where  =  ateady-atatc  drag  displacement 

Uqq  =  root  mean  square  value  c>{  the  oscillatory  drag 
displacement 

u„.  =  root  mean  square  of  the  oscillatory  lift  displace- 

OXj  . 

ment 

It  is  assumed  that  the  oscillatory  lift  and  drag  displace¬ 
ments  have  "normal"  distributions.  The  root  mean  square  values 
of  these  responses  are  equivalent  to  the  standard  deviations  (sigma) 
since  the  "means"  are  zero. 

The  use  of  3ar  (3  times  root  mean  square)  values  in  Equa¬ 
tion  46  implies  that  99.73  percent  of  the  time  the  peak  oscillating 
lift  and  drag  displacements  will  be  less  than  the  3cr  value.  Equa¬ 
tion  46  also  implies  that  the  oscillatory  components  of  the  lift  and 
drag  displacement  reach  their  maximum  at  the  same  time,  which  is 
certainly  conservative. 

5.1.2  An  Example.  An  indication  of  the  magnitude  of  bending 
moments  which  may  be  induced  by  a  strong  ground  wind  acting  on 
a  large  ballistic  missile  may  be  obtained  from  Figure  17.  In  this 
case  a  typical  intercontinental  ballistic  missile  is  assumed  to  be 
erected  on  a  surface  launcher  in  the  vertical  position  and  exposed 
to  a  uniform  steady  60-mph  wind.  The  steady- state  drag  coefficient 
is  assumed  to  be  0.55,  and  the  predominant  component  of  oscillatory 
response  in  both  the  drag  and  lift  planes  is  at  the  first  bending  mode 
frequency  of  the  missile.  The  bending  moments  at  several  points 
along  the  missile  are  shown  in  Figure  17. 

As  can  be  seen  from  these  figures,  a  strong  ground  wind  can 
easily  impose  loadings  on  a  missile  which  could  be  catastrophic  if 
not  taken  into  account  in  the  structural  design  of  the  airframe  or  if 
the  missile  is  not  otherwise  protected.  These  figures  also  indicate 
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Bending 

Moments, 

in,  lb. 

X  10' 

Distance  from  Missile  Base,  in. 

^DS 

3Mod 

3Mol 

0 

15.07 

19.54 

63.50 

72.32 

150 

9.  82 

14.  10 

44.80 

51.59 

350 

4.  84 

7.33 

23.  80 

26.75 

650 

0.  80 

0.78 

2.52 

2.97 

Figure  17.  Bending  Momenta  in  a  Typical  Missile  Exposed 
to  a  Uniform  60-mph  Wind 

the  fallacy  of  using  just  the  steady-stage  drag  response  as  a  design 
criteria,  since  the  maximum  oscillatory  components  greatly  exceed 
the  steady-state  response. 

5.2  Launching  in  a  Wind 

5.2.1  Solution  of  the  Problem.  In  the  prelaunch  condition,  the 
missile  is  bent  over  by  a  steady  wind  and  is  also  oscillating  in  its 
various  constrained  modes  because  of  the  vortex  shedding  previously 
discussed.  The  transients  that  occur  in  the  postlaunch  or  free-free 
condition  depend  very  strongly  on  the  deflected  shape  and  velocity  of 
the  missile  just  before  release.  It  would  be  desirable  to  find  the 
initial  condition  that  gives  the  highest  postlaunch  loads,  but  with  the 
missile  oscillating  in  several  modes  with  random  phase,  it  is  im¬ 
possible  to  say  a  priori  what  condition  would  give  the  largest  loads. 

It  seems  evident,  however,  that  for  any  given  constrained  mode  the 
resulting  free-free  loads  will  be  largest  if  the  missile  is  released 
when  it  is  bent  over  to  its  maximum  value.  Loads  indicative  of  the 
most  severe  postlaunch  loads  likely  to  occur  can  be  obtained  if  the 
postlaunch  transients  are  computed  using  as  initial  conditions  the 
missile  bent  into  each  of  its  constrained  modes  which  have  appreci¬ 
able  amplitude  because  of  wind-induced  oscillations.  This  means  a 
separate  computation  for  each  prelaunch  (constrained)  mode  and  a 
computation  for  the  initial  static  deflection  caused  by  steady  drag. 


-49- 


Aa  an  example  of  the  procedure  involved  in  computing  the 
postlaunch  transients,  consider  a  case  in  which  the  missile  is  bent 
into  a  deflection  curve  given  by  as  shown  in  Figure  18.  The 

function  yQ{x)  could  represent  any  of  the  displacements  discussed 
earlier,  but  for  this  example  it  will  be  assumed  It.  be  due  i.nl  /  to  ,i 
steady  wind.  The  procedure  to  be  ft.llowed  here  is  to  represent  the 
missile  in  the  postlaunch  condition  in  terms  of  its  free -free  rnodes 
0^(x),  with  generalized  coordinates  Assuming  that  the 

free-free  modes  i|,^)  have  been  obtained  from  Equation  6,  the 
problem  is  to  determine  the  generalized  coordinates  appearing  in 
Equations  9  and  10,  that  is,  solve  Equation  18  subject  to  the  ap¬ 
propriate  initial  conditions.  For  the  sake  of  convenience,  the 
pertinent  equations  are  collected  here: 

oo 

u{x,  t)  =  €^(1)  -f  xCj(t)  +  ^  ^^(x)q^(t)  (9) 

n-  1 


CD 

,1,  (x,  t)  =  C^(t)  +  ^jj(x)q^(t)  (10) 

n=  1 


EiM  =  -  M 

3x 

(4) 

<  \ 

ICAG  ‘iii.  =  V 

,  ax 

(5) 

q  +  2b  w  q  +  ui  = 

n  n 

Q 

n 

M 

n 

(18) 

The  term  involving  -  P  has  been  dropped  from  Equation  5.  The 
axial  load  P  does  not  affect  the  bending  dynamics  for  the  problem 
being  considered.  The  generalized  fc.rces  Cn  und  are  given 
by  Figure  18.  Here 
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Figure  18.  Coordinate  Syatem  of  Constrained  Missile 
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0 


(47) 


Q_ 


^  w(x)^^(x)dx  and 


and  are  independent  of  time  since  w{x)  ia  assumed  to  arise  from  a 
steady  wind.  The  solution  of  Equation  18  is 


f 


n 


_  N  -b  (jj  t  ,  . 

\e  b 

q  (t)  =  a  -  - j  ]  icos  u)  t - sin  u  t)+ 

'  I  n  ,,  2  I  1  n  u)  n  I  . .  a 

n  /Mu 

/  n  n 


I  M  u  ‘^/ 
V.  n  n  / 


(48) 


where  q  (0)  equals  a  ,  q  (0)  equals  zero,  and  b  is  assumed  to 
n ,  n  n  n 

be  small,  that  is,  the  damped  and  undamped  natural  frequencies  are 
essentially  equal.  It  only  remains  to  determine  the  constants  a^. 
To  this  end,  the  initial  bending  moment  Mq(x)  and  shear  force 
Vq{x)  caused  by  the  steady  wind  are  computed  from  (see  Figure  18) 

A. 


-Mq  = 


w(|)(|  -  x)df 


(49) 


rz 

=  w(5)dg 

Jx 


(50) 


Equations  9,  10,  4,  and  5  may  be  combined  to  yield  the  following 
expansions  for  Mq  and  Vq: 


oo 


-  Mq(x)  =  El  a^,p^-(x) 
n=  I 


(51) 


VqW  = 


n=  1 


(52) 


where  prime  denotes  differentiation  with  respect  to  x.  The  con¬ 
stants  in  these  expansions  may  now  be  determined  with  tiie  aid  of 
Equation  7,  that  is 
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'  j  M  u  m  =  n  . 

,  KACI,  ""  (53) 

X '  0  m  /  n 


Equation  51  is  now  multiplied  by  and  Equation  52  is  multiplied 

by  results  are  added  and  integrated  over  the  length  of 

the  missile  to  yield 


dx  = 


CO 

Eh|j  '  a  iij  '  dx 
,  m  /  n  n 

^  n=l 


+  )  KAG{^  '-Jj  ))  a  (0  '-Jj  )dx 
j,  "^m  ^m  Z_,  n'^^n  ^n 

^  n=  1 

(54) 

A  cr-rr.pariscr.  of  Eqi.aticns  54  and  53  yields  immediately 


a  = 
n 


M  b>  j 
n  n 


7-  ■'  '  +  V.(0  ’  -  4<  )  dx 

2  I  O  n  O'^^n  ^n  ; 


(55) 


with  the  determined,  it  is  now  possible  to  compute  the  bending 

moments  and  shear  forces  generated  as  a  result  of  launching  a 
— ii.  a  'vir.d. 


5.2.  Z  An  Example.  Because  the  moment  at  the  base  x  =  -  ij  is  the 
maximum  moment  in  the  constrained  condition,  while  the  moment  at 
this  station  is  zero  for  all  of  the  free-free  modes,  the  convergence 
of  the  initial  -i:  n.-,*n',  *  p  c-^^e.'l  in  t<- :-ms  of  free-free  modes  is  under¬ 
standably  poor.  However,  such  an  expansion  has  been  made  using  the 
bending  modes  of  an  actual  missile,  and  convergence  to  a  reasonable 
engineering  accuracy  was  obtained  using  six  bending  modes.  Figure  19 
shews  a  plot  of 

The  terms  corresponding  to  the  sloshing  coordinates  and  P  have  been 
dropped  from  Equation  7.  Experience  has  shown  that  sloshing  and  P 
are  not  important  in  the  problem  being  considered. 
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Figure  19.  Absolute  Maximum  Bending  Moments  at  Laimch 


N 

\ 


n=  1 


for  values  of  N  from  onj  to  lix.  Tht;  initial  moment  Mq(x)  is  also 
shown,  and  it  can  be  ^ee.i  tnat  conversance  is  good  for  sections  awa, 
from  the  base.  Near  the  base,  the  maximum  loads  occur  before 
launching  so  that  the  poor  convergence  here  is  not  important.  An 
indication  of  the  convergence  during  the  postlaunch  transient  is  given 
by  Figure  20,  which  shows  the  maximum  value  of  M(x,  t)  incurred 
during  the  transient  at  several  stations  of  the  missile. 

Again,  N  is  varied  from  one  to  six  in  this  figure,  and  it  can 
be  seen  that  the  additional  moment  due  to  the  sixth  mode  is  not  too 
large.  These  computations  were  made  with  no  postlaunch  damping 
(that  is,  b^  =  0).  When  the  damping  was  included,  the  convergence 
was  better,  particularly  near  the  center  of  the  missile  where  the  lower 
free-free  modes  contribute  strongly  to  the  moment.  At  these  stations, 
moments  caused  by  the  higher  :r.  des  were  damped  out  by  the  time  the 
moments  caused  by  the  lower  modes  became  in  phase.  A  comparison 
of  the  peak  postlaunch  moments  with  and  without  damping  is  shown  in 
Figure  19.  The  amount  of  damping  used  in  the  analysis  was 


b  =0.02 


Figure  19  also  indicates  the  importance  of  the  postlaunch 
transient.  Notice  that  in  the  prelaunch  co.ndition  the  bending 
moments  near  the  tip  of  the  missile  are  quite  small,  whereas  during 
the  postlaunch  transient  the  bending  moments  become  quite  appreci¬ 
able  because  of  the  unloading  wave  induced  by  the  sudden  release. 

The  characteristi  dip  in  the  moment  curve  at  station  0.44  of  Figure  19 
is  caused  by  a  sudden  reduction  in  the  flexural  rigidity  of  the  missile 
at  this  station.  These  postlaunch  moments  designed  the  forward  sec¬ 
tion  of  the  missile  of  this  example.  Details  of  this  analysis  can  be 
found  in  Reference  22. 
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56 


Absolute  Sum  of  Peak  Mode  Responses  of  Bending 
Moments  Due  to  Sudden  Release  from  Initial 


5.3  Effocta  of  Thrust  Dulld-up 


In  this  section  the  problem  of  determining  the  dynamic  loads 
generated  by  an  engine  start  transient  is  considered.  This  problem 
differs  from  those  discussed  in  previous  sections  in  that  the 
longitudinal  rather  than  the  bending  dynamics  of  the  airframe  is  in¬ 
volved. 

5.3.  1  Theoretical  Considerations.  Consider  a  missile  erected  on 
and  attached  to  a  launch  stand  (see  Figure  21).  At  time  t  =  0  the 
engines  are  started,  generating  a  time  dependent  thrust  F(t).  The 
problem  is  to  determine  the  response  of  the  airframe  to  F(t). 

Experimental  data  have  confirmed  that  the  longitudinal  dy¬ 
namics  of  a  missile  airframe  are  adequately  represented  by  simple, 
one-dimensional,  beam  theory.  Therefore,  the  longitudinal  respons 
of  the  missile  is  determined  from 


a^{x)  =  9^u(x,  t) 

ax^  at^ 


subject  to  the  boundary  conditions.  Here 


8u(o.  t)  _  g  au(S,  t) 
ax  '  ax 


-  SKTe  [^<*>  *  •'] 


where 

a^(x) 

E 

A(x) 

m{x) 
u{x,  t) 

X 


EA(x) 

m<x) 

modulus  of  elasticity 

cross-sectional  area  of  the  load  carrying 
portion  of  the  airframe 

mass  per  unit  length  of  the  airframe 

axial  displacement  of  the  airframe 

axial  coordinate  measured  from  the  missile  tip 

spring  constant  of  the  launcher 


e  and  Mass-Spring  Analog 


A  direct  integration  of  Equation  56  is  usually  quite  difficult, 
especially  when  the  effects  of  liquid  propellants  are  included. 
Therefore,  rather  than  attempt  a  direct  solution  of  Equation  56, 

It  is  far  more  convenient  to  replace  the  continuous  airframe  by 
a  lum,jcd  j>ring-masB  analog.  A  typical  decomposition  of  an  air¬ 
frame  into  a  spring-mass  system  is  shown  in  Figure  21. 

The  general  procedure  used  to  arrive  at  such  an  analog  is 
to  divide  the  airframe  into  a  number  of  sections,  replacing  the 
distributed  mass  of  each  section  by  a  single  mass  placed  at  the 
center  of  gravity  of  each  section;  the  masses  are  joined  by  springs 
whose  stiffnesses  are  equal  to  the  stiffness  of  the  airframe  struc¬ 
ture  between  the  masses.  The  inclusion  of  propellant  effects  re¬ 
quires  special  treatm.ent.  The  hydrostatic  pressure  created  by 
the  effective  weight  of  accelerating  propellants  causes  the  propel¬ 
lant  tanks  to  bulge,  which  in  turn  shifts  the  center  of  gravity  of  the 
propellant  mass.  A  computation  of  these  events  leads  to  expressions 
for  effective  spring  constants  of  the  form  given  by  K^,  and  Kj 

of  Figure  21.  The  expressions  for  these  constants  are  derived  in 
Section  1  of  the  Appendix. 

The  equations  of  motion  of  the  spring-mass  analog  can  be 
written  by  inspection.  Some  of  these  equations  are 


M  X  +  K  (x 
n  n  n  n 


•  1  n  n 


■  ^n-l> 


0 


M 


,x  ,  +  K.,{x  ,  -  X  )  +  K^{x  ,  -  X.)  +  K  (x  ,  -  X  ) 

n-1  n-1  7  n-I  p  5  n-1  4  n  n-l  n 


+  C  (i  ,  -  X  )  =  0 
n  n-1  n 


M  X  +  K,(x  -  X.)  +  K_(x  -  X  ,)  =  0 

p  p  6'  p  4  7'  p  n-1 
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M 


^  ^6^^4  *  ’‘p^ 

+  VP  =  ° 

^x^  +  Kjlx^  -  Xj)  t  Cjix^  -  Xj)  +  ^^{x^  -  x^)  +  ^^(><2  '  x^) 


+  Vz  "  ° 


-e  W 

"2-1^ 


^2^2  ^  ^1^^2  "  ^  ^  ~  ^4^ 


+  K^(x2  -  x^)  +  W  -  0 


W 


where 


n 


W  =  )  rn.g 
i-^1 


The  limits  on  the  last  two  equations  of  58  account  for  the 
possibility  of  the  missile  separating  from  the  launcher.  The 
solution  of  these  equations  provides  axial  loads  F^,  for  example, 


F 

n 


K  (x  -  X 
n  n  n- 


1  n  n 


and  accelerations  x  for  different  sections  of  the  missile.  Since 
n 

the  coordinates  originate  from  the  equilibrium  position,  the  absolute 
acceleration  of  any  mass  will  be  +  g,  where  g  is  the  accelera¬ 
tion  caused  by  gravity.  Similarly,  the  total  axial  load  at  station  n 
equals  K(x  -x  ,)-mg  +  C(x  -x  ,),  This  choice  of  co- 
ordinates  is  made  in  order  to  eliminate  the  necessity  of  computing 
intial  displacements  and  setting  them  as  initial  conditions  on  a 
computer.  The  initial  displacements  for  the  equations  here  are  all 
zero,  whereas  if  all  spring  deflections  were  measured  from  their 


(58) 
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froQ  lonfjth,  the  deflection  of  each  spring  under  the  weight  of  the 
masses  above  it  would  have  to  be  used  to  determine  initial  condi¬ 
tions.  In  a  more  complicated  model  with  springs  spanning  across 
several  masses,  solution  of  the  static  problem  is  more  tedious  so 
that  the  use  of  coordinates  measured  from  equiUbrium  is  even 
more  strongly  indicated. 


In  many  engine  thrust  transient  response  problems  the  effect 
of  damping  in  the  missile  is  very  important.  Hence,  an  estimate 
of  the  modal  damping  {possibly  an  experimental  value)  for  the 
structure  is  necessary  to  determine  the  values  of  dashpots  C^, 
which  are  located  between  masses  as  shown  in  Figure  21.  If  the 
modal  damping  b^  in  the  missile  is  small,  it  can  be  shown 
(Re  ference  15)  that  a  good  approximatior  for  b^  is  given  by 


2u>  M 
n  n 


where 


b  =  percent  of  critical  viscous  damping 
in  nth  mode 


Cj  =  dashpot  constant  at  station  j 
P  =  total  number  of  dashpots 


relative  displacement  in  dashpot  j, 
for  the  nth  mode 


=  circular  frequency  of  nth  mode 

M  =  generalized  missile  mass  for  the  nth 
mode 

N  2 

iTi 


N  =  total  number  of  masses 


(59) 
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Equation  59  provides  the  necessary  values  of  the  dashpot  constants 

CO  maintain  the  d-:  sired  modal  damping.  A  set  of  P  simultaneous 

equations  are  formed  with  P  dashpot  constants,  C.  as  unknowns. 

J 

Equation  58  can  be  written  in  matrix  form  as  follows; 

[aJ{x>  +  [bJ'x)  +  [cj/x)  =  0  (60) 

mass  matrix 
damping  matrix 
stiffness  matrix 
a  column  matrix 

For  purposes  of  computing  longitudinal  modes  and  natural  fre¬ 
quencies,  the  damping  is  assumed  to  be  negligible,  and  the 
time  dependence  of  the  displacements  is  sinusoidal  with  fre¬ 
quency  «.  Equation  60  then  takes  the  form 


fAl  = 
[b]  = 

s  . 


([c]  -  u,2[a])(x)  =  0  (61) 

Equation  61  is  the  classical  form  for  an  eigenvalue  problem  in 
matrix  notation.  The  solution  of  61,  yields  the  natural  frequencies 
and  mode  shapes  p^{x).  With  the  mode  shapes  available,  the 
response  of  the  spring -mass  analog  can  be  written  as 


N 

u{x,  t)  =  ^  ^>^(x)q^(t) 

n=  1 


with  the  generalized  coordinates  q  (t)  given  by 


q  (t)  +  2b  u)  q  (t)  +  (aj  q  (t)  = 


F{t)^^(xJ 
M 


(62) 


(63) 


■Equations  62  and  63  are  valid,  independent  of  whether  the  missile  is 
resting  on  or  has  left  the  launcher,  the  only  difference  being  that  the 
mode  shapes  are  different  in  each  case. 

Equations  62  and  63  may  not  represent  the  most  convenient 
form  of  the  solution.  In  fact,  it  is  often  much  faster  to  solve 
Equation  58  directly  on  an  analog  computer  rather  than  attempt  a 
modal  solution. 

5.3.2  An  Example.  Figure  22  shows  the  effects  of  damping  and 
release  time  on  missile  axial  loads  caused  by  launch  and  engine 
thrust  transients.  The  axial  loads  occur  in  a  large  ballistic  mis¬ 
sile  just  after  launch  when  the  release  time  is  specified  at  a  per¬ 
cent  of  the  full  engine  thrust.  A  typical  engine  thrust  build-up 
cur'/e  is  shown  in  Figure  23  with  release  times  of  the  missile 
which  correspond  to  the  cases  of  Figure  22.  To  compare  cases 
B  and  D,  Ll\e  differences  in  the  loads  are  caused  by  the  transients 
of  the  engine  thrust  which  are  generated  in  part  by  the  rate  of 
build-up.  In  case  D  these  transients  are  nearly  damped  out,  and 
hence  for  this  case  the  results  may  be  obtained  by  simply  applying 
a  step  load  to  the  end  of  a  free-free  missile  of  magnitude  T  -  W, 
where  T  is  the  100  percent  thrust  value  and  W  is  the  total  weight 
of  the  missile.  To  compare  cases  A  and  B,  the  differences  in  the 
axial  loads  is  mainly  caused  by  the  differences  in  the  rigid-body 
accelerations.  The  effect  of  damping  is  clearly  indicated  in  cases 
B  and  C.  Although  the  axial  loads  of  Figure  22  indicate  that  it  is 
advantageous  to  release  the  missile  at  the  lowest  thrust  level,  the 
combination  of  axial  and  lateral  loads  may  indicate  the  necessity 
for  a  different  time  of  release.  An  over-all  loads  picture  must 
be  examined  before  a  release  rime  can  be  specified.  The  lateral 
loads  will  dominate  the  axial  loads  in  certain  sections  of  the  mis¬ 
sile  and  vice  versa.  Therefore,  if  the  lateral  loads  are  predominant, 
which  is  usually  the  case  at  sections  near  the  tip  of  the  missile,  it 
may  be  necessary  to  launch  at  a  time  when  some  of  the  contributors 
to  lateral  loads  are  damped  out,  that  is,  effects  of  engine  thrust 
differential,  which  have  not  been  treated  here. 
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Case 

Missile  Station 

Maximum  Compressive 

Axial  Load 

A: 

125  tip 

1.0^ 

Release  at  77%  of  full  thrust 

250 

4.  5 

=  1% 

400 

12.2 

n 

800 

38.  8 

n  =  1 ,  2,  3,  4 

B: 

125 

2.  1 

Release  at  100%  of  thrust 

250 

11.1 

b  =  1% 

400 

27.  5 

n 

800 

64.  5 

C: 

125 

3,8’ 

Release  at  100% 

250 

18.0 

b  =  0 

400 

45.0 

n 

800 

73.  5 

D: 

125 

1.9 

itc lease  after  a  1- second  hold 

250 

9.7 

b^  =  1% 

400 

24.  0 

n 

800 

59.0 

^Loads  are  normalized  to  {b)  and  occur  after  release  (postlaunch). 


Figure  E2.  Effect  of  Damping  and  Release  Time  on  Missile  Axial 
Loads  due  to  Launch  and  Engine  Thrust  Transients 


Figure  23.  Engine  Thrust  Build-up  Curve 
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5.4  Mlfifllla  Response  duo  to  Ground  Shock 

Many  dynamic  problems  are  encountered  before  the  missile 
is  airborne,  for  example,  in  transport  and  ground  handling,  how¬ 
ever,  this  section  will  be  limited  to  the  discussion  of  missile 
dynamic  response  due  to  ground  shock. 

There  are  many  complex  problems  encountered  in  determin¬ 
ing  the  responses  of  missiles  due  to  ground  shock  and  in  determin¬ 
ing  the  characteristics  of  shock  Isolation  systems  which  attenuate 
these  responses.  Some  investigators  have  examined  these  problems 
by  representing  a  continuous  structure  by  a  single  degree  of  free¬ 
dom  system.  This  system,  however,  may  be  inadequate  or  the 
responses  such  as  displacements,  moments  and  shears  may  be 
misleading.  Therefore,  the  more  difficult  problems  associated 
with  a  continuous  structure  should  oe  investigated.  A  method  of 
analysis  for  these  prc  b'  .  using  a  "shock  spectra"  technique 
is  given  ‘n  Refc^i-nr.*  ^  •  lud 

In  this  section  two  aspects  ■>!  thci'a  p’^oblems  are  considered. 
They  are: 

a)  The  specification  of  the  shock  environment 

b)  The  dynamic  response  of  the  missile  to  the 
environment. 

The  ground  shock  environment  may  be  specified  in  various 
ways,  such  as  by  soil  pressure,  displacement,  permanent  set, 
velocity  or  acceleration.  All  of  these  may  be  of  use  in  particular 
applications.  For  the  purpose  of  this  section  attention  is  given  to 
the  displacement  shock  spectra  (Reference  25). 
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Although  thure  are  certain  practical  difficulties,  It  is 
theoretically  possible  to  determine  the  response  of  a  structure 
to  shock  if  the  acceleration  time  input  and  the  dynamic  charac¬ 
teristics  of  the  structure  in  terms  of  the  principal  modes  fre¬ 
quencies,  and  damping  are  known.  It  is  also  possible  to  estimate 
upper  bounds  of  response  by  means  of  the  shock  spectra  and  the 
dynamic  characteristics  of  the  structure.  This  latter  procedure 
will  be  discussed  in  the  determination  of  maximum  stresses  and 
displacements,  i*.  is  p'^inied  out  to  tLe  reader  that  this  method 
of  determining  the  upper  b  ..u.d  response  of  a  structure  due  to 
transient  excitation  is  aiso  i  very  useful  tool  for  problems  other 
than  ground  shock.  It  has  important  applications  ir.  the  specifica¬ 
tion  of  the  type  and  intensities  of  shocks  a  structure  must  resist,' 
and  in  the  estimation  of  tiie  strength  of  a  structure  in  resistance 
to  a  variety  of  shocks.  Such  problems  arise  in  the  packaging  of 
electronic  equipments  for  transportation,  in  the  resistance  of  a 
structure  to  earthquake,  in  the  landing  of  an  aircraft,  in  the 
starting  shock  of  a  rocket,  and  in  a  large  variety  of  impact  loads. 

To  discuss  the  response  of  missiles  to  ground  shock,  the 
missile  is  considered  in  a  vertical  position  attached  to  a  launch 
stand  which  is  protected  from  air  blast  but  is  subjected  to  ground 
shock.  (See  Figure  24.)  The  ground  shock  in  terms  of  the  soil 
acceleration  as  a  function  of  time  consists  of  two  major  compo¬ 
nents:  a  vertical  component  and  a  component  directed  radially 
from  ground  zero  where  the  blast  initiates. 

The  dynamic  characteristics  of  the  structure  are  con¬ 
veniently  expressed  in  terms  of  the  principal  modes  and  fre- 
cuenci  rs,  and  Etruf:'.ural  damping.  For  the  missile-stand 
combination  the  transverse  modes  are  very  similar  to  those  .of 
an  elastically  supported  cantilever  beam.  The  displacement  of 
any  point  of  the  structure  can  then  be  taken  as  the  superposition 
of  the  mode  contributions. 
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Figure  24,  Mis alle  — Launcher 
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The  diaplacomant  of  the  miaaile  relative  to  the  ground,  as 
shown  in  Figure  24  is  given  by, 


u(x,  t) 


(64) 


where 

u  =  ..iaj;.la.-eii. en.  of  the  structure  at  point  x 
relativ.  to  the  ground 

=  generalized  coordinate  of  the  mode 
=  value  of  the  mode  at  position  x 


For  an  acceleration  input,  the  equation  of  motion  for  the  n^^  mode 
is. 


q_+2twq  +aj  q„  = 
m  ’n  nm  n  ^n 


-  a{t) p^jjdx 


dx 


(65) 


OJ 

n 


P 


a(t) 


ratio  of  damping  to  critical  viscous  damping 
in  the  nth  mode 

circular  frequency  of  the  n^^  mode 

mass  distribution  per  unit  length 
acceleration  of  the  ground 


The  right  hand  side  of  this  equation  represents  the  generalized 

f  **  2 

force,  a(t)  /  dx,  divided  by  the  generalized  mass,  /  p^  dx. 

Ji  ^  r  /r  2  Jjt 

.For  convenience  let  y  =  /,  p0  dx/  L  pd  dx.  This  term  has  been 

n  Jt  / Jt 

called  the  kinematic  or  distribution  factor  (Reference  24). 
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The  solution  to  the  foregoing  equation  of  motion  for  small  damping 

and  for  initial  conditions  of  a,  q  ,  and  q  equal  to  zero  at  t  =  0, 

n  n  * 

is. 


w  (t-T) 

q  (t)= - I  a(T)e  sin  u  (t  -  T)dT 

“n  -JO  n'  ' 

which  when  substituted  into  Equation  64  determines  the  response 
of  the  structure. 

There  are  several  practical  difficulties  in  determining  the 
response  of  the  structure  in  the  foregoing  manner.  Some  of 
these  difficulties  come  from  the  evaluation  of  the  integral  with 
practical  acceleration  records  a(t);  others  arise  from  the 
characteristics  of  the  acceleration  records  such  as  statistical 
scatter,  damping,  and  determination  of  low-frequency  components 
These  difficulties  may  be  present  in  problems  other  than  ground 
shock. 

It  is  sometimes  convenient  to  think  of  the  relationship  of 
the  generalized  coordinates  to  quantities  that  can  be  measured 
directly  without  a  knowledge  of  the  ground  accelerations.  In 
effect,  the  structural  system  for  which  dynamic  responses  are 
to  be  determined  is  replaced  by  a  series  of  single-degree-of- 
freedom  systems.  For  example,  the  analog  for  the  motion  of 
each  mode  of  the  missile  on  the  stand  can  be  represented  by  a 
cantilever  beam  tuned  to  the  principal  frequency  for  the  mode 
of  the  missile.  Hence,  assuming  an  idealized  single-degree-of- 
freedom  system  such  as  a  cantilever  with  an  end  mass,  the 
equation  of  motion  for  the  mass  is. 


Thia  equation  has  the  aame  form  aa  Equation  (65),  ao  that  if  the  motion 
o£  the  maas  of  the  alngla-degree-of-freedom  ayatem  la  known  (by 
meaauremont  or  other  moana),  the  generalized  coordinate  for  the  nth 
mode  of  the  structure  la  given  by 

q  =  Y  Q 
^n  n  n 

and  the  structural  response  is 


(66) 


Now  the  motion,  Q(t),  of  the  mass  ia  often  difficult  to  measure 
as  a  complete  displacement-time  history.  Fortunately,  it  is  possible 
to  make  some  predictions  about  the  maximum  value  of  the  response  of 
the  structure  from  a  knowledge  of  only  the  peak  values  of  the  response 
of  the  3 ingle -degree -of-freedom  system,  which  are  relatively  easy  to 
calculate  if  the  pulse  input  is  known  or  to  measure  by  means  of  reed 
shock  gages  such  as  in  the  case  of  ground  shock  input.  (See  Reference  24.) 

A  shock  gage  may  consist  of  a  set  of  single-degree-of-freedom 
systems,  as  shown  in  Figure  25-  The  frequency  spectrum  of  the  peak 
displacements  of  the  masses  relative  to  the  base  which  is  being  ac¬ 
celerated  is  called  the  displacement  shock  spectrum  and  is  given  by 


D(u)  = 


ax 


t  o 


^  ^  a(  T)e  sin  w  (t  -  t)  dr 


Schematically,  D  and  Q  for  a  frequency  u  are  shown  in  Figure  25  with 
the  associated  ground  displacement  curve.  A  typical  horizontal  dis¬ 
placement  shock  spectrum  is  shown  in  Figure  26. 
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Figure  26.  HorlzontaJ.  Displacement  Shock  Spectrum 


The  greatest  value  oC  the  response  in  Equation  (66)  would  occur  if 
all  the  peaks  of  the  components  occurred  simultaneously  and  in  the  same 
phase.  Since  the  peak  value  of  Q(t)  is  D,  then 


u{x) 


Y  D(cj  )<p  (x) 
'n  '  n'^^n' 


Having  an  expression  for  the  displacement,  other  quantities  of  interest, 
such  as  bending  moment,  can  be  obtained.  The  expression  for  bending 
moment  is. 


Y  D(u  )m  (x) 
'n  n  n 


where  m  (x)  is  the  bending  moment  at  position  x  in  the  nth  mode, 
n 

These  estimates  of  response  are  upper  bounds;  however,  it  has  been 
shown  by  several  investigators  (References  23,  25,  26  and  27)  that  for 
shock  inputs  such  as  earthquake  or  ground  shock,  the  estimate  is  not 
more  than  15  percent  conservative.  The  conservativeness  of  the  re¬ 
sponse  when  using  single  pulse  shock  inputs  is  discussed  in  detail  in 
Reference  25. 


Such  a  "shock- spectra"  technique  of  analysis  is  also  a  very 
valuable  tool  in  determining  the  response  of  a  missile  during  engine 
thrust  build-up  if  the  boundary  conditions  do  not  change.  For  example, 
if  a  missile  is  constrained  in  a  launcher  or  if  it  is  free-free  as  in  the 
case  of  a  space  vehicle  at  the  time  of  engine  ignition  after  a  coast 
period,  an  upper  bound  of  response  may  be  obtained  by  using  a  modal 
analysis  if  the  input  force  (engine  thrust)  is  represented  by  its  "shock 
spectrum" . 
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APPENDIX 


MODELS  FOR  THE  LONGITUDINAL  VIBRATIONS  OF  LIQUID 
PROPELLANT  MISSILES 


Equations  for  a  Simple  Cylindrical  Tank 


A  lumped  parameter  model  is  constructed  to  approximate  the  first 

mode  of  vibration  of  a  liquid-filled  elastic  tank  which  is  considered  to  be 

a  tandenl  part  of  a  missile  as  in  Figure  21.  The  propellant  tank  is  shown 

in  Figure  27.  The  force  F  In  the  figure  is  the  force  in  the  tank  walls  due 

to  the  action  of  masses  above  the  tank  under  consideration.  is  the 

ullage  pressure  in  the  tank.  Other  quantities  are  evident  from  the  figure. 

The  axial  stress,  cr  ,  and  radial  stress,  <r  ,  in  the  tank  wall  are; 

X  y 

-F  ^o  ^ 

''  2^?RT  (A.1) 


riuio  AccrwenAnoH  a 


Figure  27.  Propellant  Tank 
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(A.  6) 


w=  -R,y 


P  +  pa  (h  -  x)  + 


vF 

2ttR^. 


In  solving  a  longitudinal  transient  problem,  the  stresses  due  to 
the  static  ullage  pressure,  P^,  can  be  added  separately  to  the  dynamic 
stresses.  Dropping  these  static  terms,  the  displacement  of  the  section 
of  the  tank  at  the  free  fluid  surface  relative  to  the  base  is. 


u(h)  =  - 


R  r  Fh 
IzirR^ 


+ 


(A.  7) 


The  total  mass  of  the  propellant  and  a  lumped  spring  constant  for  the 
tank  are  respectively. 


M  =  itR  hp 
P 

AE  _  ZirRtE 
■  h  "  h 


The  displacement  of  Equation  {A..  7)  can  be  written  as, 

vM 


u(h)  = 


F  ■  P 
K  ■  K 


X 

o 


where; 


(A.  8) 
(A.  9) 


(A.  10) 


Xj  =  coordinate  of  the  tank  position  at  free  surface  of 
the  liquid  measured  from  static  equilibrium. 


X  =  coordinate  of  the  tank  base,  measured  from  static 
equilibrium 


Assume  now  that  the  fluid  can  bo  treated  as  a  rigid  body,  and 
that  its  displacement  Is  the  vertical  displacement  of  tlie  center  of  gravity 
of  the  fluid  which  occurs  when  the  tank  walls  expand  from  w  =  0  to 
w  =  w(x)  as  given  by  Equation  (A.  6).  Neglecting  second  order  volumes, 
this  distance  (see  Figure  28)  is, 

1 

X  -  X  = - r  2irRw(x)  dx 

P  “  irR^-Jo 

3M  „„ 

=  -~K 

where,  x^  =  coordinate  of  fluid  eg,  measured  from  static  equilibrium 

Eliminating  F  from  Equations  (A.  10)  and  (A.  11)  and  noting  that 
a  =  Xp,  one  obtains; 

3M  r  vM  • 

X  -  X  =  -  P-  X  +  V  (x ,  -  X  )  +  ■  X  (A.  1 2) 

p  o  ZK  p  [  1  o  K  p  \  I 


This  equation  can  be  rewritten  using, 


Cx,  -  x)sCx,  -  x)-{x  -x) 

1  o  Ip  p  o 


,  (2  -  2v)  K  ,  ,  2vK  ,  ,  _  „ 

M  X  +  - y  (x  -  X  ) - y  (x,  -  X  )  -  0 

P  P  (3  -  2V  )  P  °  (3  -  2v  )  ^  P 


(A,  13) 


Notice  that  Equation  (A,  13)  corresponds  to  Figure  29  if  one  takes. 


(A,  14) 
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Figure  28.  Displacement  of  Fluid  Figure  29.  Longitudinal  Model 
due  to  Tank  Wali  of  Propellant  Tank 

Expansion 


The  spring  has  been  included,  because  neither  Equation  CA.  10)  nor 
(A.  11)  gives  the  same  deflections  as  found  in  the  figure  if  =  0.  To 
find  K^,  set  =  0  and  make  the  figure  match  the  equation, 

■  ’'o  =  ■  ^  (A,  15) 


Using  Equation  (A.  14)  the  force  in  spring  is, 


(A.  16) 


Also,  the  deflection  of  spring  K-,  using  Equation  (A.  10  with  x  =  0,  is 

J  P 


X,  -  X 

1  o 


F 

K 


Equilibrium  of  the  lower  bar  in  Figure  29  under  thoae  forces  requires 
that, 


(A.  17) 


CA.  18) 


In  Figure  21,  therefore,  by  using  Equations  (A.  14)  and  (A.  18) 


B.  Model  for  a  Cylindrical  Tank  With  a  Flexible  Lower  Bulkhead 

If,  instead  of  having  a  rigid  base,  as  implied  for  the  tank  in 
Section  A,  the  tank  lower  bulkhead  has  ah  equivalent  spring  constant 

the  model  in  Figure  29  must  be  modified  as  shown  in  Figure  30. 
It  is  evident  that  the  static  deflection  of  the  liquid  mass  relative  to  the 
reference  plane  is  merely  that  due  to  the  flexible  tank  plus  that  due 
to  the  flexible  lower  bulkhead,  so  that  the  spring  from  the  mass  to  the 
reference  plane  is  simply  these  two  effective  springs  in  series,  as 
shown.  When  a  force  F  is  applied  to  the  tank  skin  from  above,  the 


()*  flrc'lii)n.M  of  anrl  with  rcHpcct  to  tho  ruferenco  plane  arc  the  same 
aa  if  tlic  hiilkhoad  were  rigid  because  no  further  load  l:i  carried  by  the 
bulkhead  (i.e.  it  carries  only  the  fluid).  Thus, 


(A.  19) 


as  before,  the  new  spring  constants,  and  are  now  determined  so 
that  these  deflections  indeed  result.  From  Figure  30  and  Equation  (A.  19) 
one  can  write. 


from  which 


Kj  =  K  -  i.K^ 


(*r.  20) 


Figure  30.  Longitudinal  Model  of  Propellant  Tank  With 
Flexible  Lower  Bulkhead 
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From  Figure  30,  is  computed  to  be, 


/“  2  -  2 

(3  -  Zv^J 

(3  -  Zv^J 

Also  from  Figure  30  and  Equation  (A.  19)» 


BO  that. 


(A.  21) 


C,  Model  for  a  Semi-Monocoque  Cylinder 

If  the  cylinder  is  constructed  from  a  skin- stringer  design  as  shown 
In  Figure  31  it  ia  assumed  that  the  stringers  have  negligible  radial 
displacement  but  that  the  skin  will  displace  both  radially  and  longitudinally. 
The  L'.ii.il  force  T.'th  is  either  through  the  skin  or  the  stringers,  therefore 
the  stringers  act  in  parallel  with  the  skin.  The  model  for  the  skin  ia 
similar  to  the  case  in  Section  A.  ,  but  now, 


K 

m 


h 


and  the  stringer  spring  constant  is. 


K 

B 


(AE)  . 

_ stringer 

h 
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Figure  31.  Section  of  Skin- Stringer  Cylinder 
and  Longitudinal  Model 


( 


D,  Model  for  a  Sem£-Monocoque  Cylinder  With  Partially  Buckled  Skin 

A  longitudinal  model  of  a  skin-stringer  propellant  tank  with  a 
partially  buckled  skm  is  necessary  for  some  of  the  ballistic  missiles. 
As  in  the  previous  section  it  is  assumed  that  the  axial  load  is  trans¬ 
mitted  by  the  stringers,  K^,  and  by  the  effective  (unbuckled)  skin, 
where  the  latter  is. 


(AE) 


skin 


effective 


The  spring-mass  model  to  represent  this  case  is  shown  in  Figure  32. 
♦  ^ 

The  constants  and  K  are  such  that, 


X,  -  X  = 

1  o 


^MB  +  ^s 


(A.  22) 
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and  tljo  propellant  maBS, 
carried  by 

K^.  l.e.  . 


deflects  if  a  portion  of  the  load  F,  say 
is  transmitted  by  the  unbuckled  skin, 


X 

P 


X 

o 


vF 


MB 


-K 


M 


(see  Equation  A,  15)  but. 


therefore, 


^MB 

^MB 


^MB  + 


where  Kj^ 


(AE) 


skin 


X  -  X 

P  o 


,F  K 


MB 


MB 


Equations  (A.  22)  and  (A.  23)  are  the  constraints  to  be  satisfied. 


Assuming,  as  before,  that  x  =0  then  from  Figure  32  and 
Equation  (A.  22), 


therefore, 


(A.  23) 


(A.  24) 
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Figure  32.  Longitudinal  Model  of  Semi-Monncoque 
Cylinder  with  Partially  Buckled  Skin 


From  Figure  32  and  Equation  (A.  23), 

r  ^  -  Iv 


h  ■  -o] 


13  -  2t;' 

“  7  ^ 

V 


**SiB 


therefore. 


-.■iriT) 


K 


MB 


Substituting  Equation  (A.  25)  into  Equation  (A.  24)  gives. 


K--fc5- 

"  K.,  -V 


^M  ^MB 


^M  ■  ''^MB 


(A.  25) 


(A.  26) 


The  foregoing  derived  spring  constants  for  Figure  32  are  shown  in 
Figure  33. 
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